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We investigate if the hemispherical asymmetry in the CMB is produced from “asymmetric” excited 
initial condition. We show that in the limit where the deviations from the Bunch-Davies vacuum is 
large and the scale of new physics is maximally separated from the inflationary Hubble parameter, 
the primordial power spectrum is modulated only by position dependent dipole and quadrupole 
terms. Requiring the dipole contribution in the power spectrum to account for the observed power 
asymmetry, A = 0.07 =b 0.022, we show that the amount of quadrupole terms is roughly equal to 
A 2 . The mean local bispectrum, which gets enhanced for the excited initial state, is within the 
la bound of Planck 2015 results for a large field model, /nl — 4.17, but is reachable by future 
CMB experiments. The amplitude of the local non-gaussianity modulates around this mean value, 
depending on the angle that the correlated patches on the 2d CMB surface make with the preferred 
direction. The amount of variation minimizes for the configuration in which the short and long 
wavelengths modes are around the preferred pole and Ifel ~ |fe~io| \ki\ ~ |fe| ~ |Zfe 25 oo| with 
fmJ 1 ~ 3.64. The maximum occurs when these modes are at the antipode of the preferred pole, 
/nl x ~ 4.81 . The difference of non-gaussianity between these two configurations is as large as 
~ 1.17 which can be used to distinguish this scenario from other scenarios that try to explain the 
observed hemispherical asymmetry. 


I. INTRODUCTION 

Inflation, despite successfully explaining the general 
pattern of observed cosmic microwave background radi¬ 
ation (CMB) \l]. fails to explain few anomalies at large 
scales. Some of these anomalies, which were previously 
observed in the WMAP data [2], persist even in the latest 
Planck data, even though their statistical significances 
might not be substantial 0. Such anomalies in general 
break the statistical isotropy of the CMB and can be 
modeled phenomenologically as 

AT(i) = AT iso (x)(l + M(x )), (1) 

where T iso is the isotropic part of the temperature fluc¬ 
tuations, where x is where you look in the sky. In 
particular, there seems to be hemispherical asymmetry 
consistent with the existence of a dipolar modulation 
term, M(x) = Ax • n in the Planck data with ampli¬ 
tude A « 6 — 7% on scales 2 < Z < 64 [IJ 09 0]. As 

stated before, x is where you look in the sky and n is 
the preferred direction. The asymmetry seems to fade 
away at smaller scales, especially for l > 600 [5j |6]. The 
asymmetry is more than twice as large as the expected 
asymmetry due to cosmic variance, A « 2.9%. Similar 
asymmetry can arise from a dipolar term in the primor¬ 
dial power spectrum [7: or from a phenomenological x- 
dependent modulation of the primordial spectrum [B 

Vs = ?iso (1 + 2A(x • n)) . (2) 

Various proposals have been offered to explain this asym¬ 
metry. Some considered the effect of long wave-length 
super-horizon mode on the sub-horizon power spectrum 


012! through Grishchuk-Zekdovick effect [10|, which re¬ 
quires non-negligible local non-gaussianity to correlate 
the long and short wavelength modes. Also noncommu- 
tative physics at Planck scale [7], isocurvature perturba¬ 
tion mi, non- Gaussianity PUCE], domain walls m or 
running of scalar spectral index [15] have been suggested 
as mechanisms explaining the observed asymmetry in the 
power spectra. In principle, higher order multipoles can 
also contribute to 0, and respectively §. 

The main goal of this paper is to design a scenario in 
which the observed dipole asymmetry is realized, assum¬ 
ing that the initial condition for fluctuations has a small 
anisotropic position-dependent asymmetric partQ Such 
asymmetric-contribution in the initial condition could be 
the effect of preinflationary patch which was probably 
highly inhomogeneous and anisotropic or the effect of 
parity violating terms in the fundamental theory higher 
than the energy scale of inflation, which was also position 
dependent within the inflationary patch. The quantum 
vacuum state of the universe is thus not assumed to ini¬ 
tially be aligned with the Bunch-Davies vacuum which is 
the standard choice. In reference m non-Bunch Davies 
vacuum was also considered as a possible reason for the 
power asymmetry, but their mechanism, that was based 
on coupling of modes in an isotropic vacuum, was differ¬ 
ent. 


1 By “asymmetric”, we mean that the initial condition for the 
scalar perturbations is not invariant under the the transforma¬ 

tion x —> —x within the horizon patch (by horizon patch we 
mean the patch that becomes the size of the current observable 
universe after inflation and subbsequent stages of cosmological 
evolution.). 
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Fixing the amount of the asymmetry from the ob¬ 
served hemispherical asymmetry, we also predict a non- 
negligible quadrupole contribution to the primordial 
power spectrum 

^5 = J’iso (l + 2A(x ■ h) + B(x ■ h) 2 ) . (3) 

with B within the interval 

0.0025 <B ~ A 2 < 0.008. (4) 

In our scenario, higher order multipole contributions to 
the primordial power spectrum will not only be sup¬ 
pressed by higher powers of A , which is small, but also 
by negative powers of Nk 1, where Nk is the number 
of quanta in the initial excited state which is related to 
the second Bogoliubov coefficient /3k through the rela¬ 
tion Nk = \/3k \ 2 . As shown in [18], one can start from ex¬ 
cited initial states with large occupation number, Nk 1 
without violating the bounds on backreaction. 

The structure of the paper is as follows. First we re¬ 
view the formalism of excited initial condition, show¬ 
ing how dipole and quadrupole terms could be gener¬ 
ated from a asymmetric excited initial condition. As ex¬ 
pected, the local configuration is enhanced for such ex¬ 
cited initial states with an amplitude which is within the 
la bound of Planck data. However, in addition, one finds 
an angular-dependent modulation that depends on the 
direction that each mode is located at and the angle that 
it makes with the preferred direction. The amplitude of 
modulation minimizes for the local configuration around 
the preferred pole, in which the short wavelength modes 
are the smallest one that could be probed by Planck, 
l ~ 2500 and the longest one corresponds to l ~ 10 at 
which the cosmological variance is small. For the same 
local configuration which is at the antipode pole the local 
non-gaussianity reaches its maximum. We conclude the 
paper at the end. 


II. ASYMMETRIC EXCITED INITIAL STATES 

As quite well-known, the predictions of inflationary 
models for the CMB spectrum depends on the initial 
state of the quantum perturbations as well as the spe¬ 
cific details of the model. The standard lore is that these 
perturbations embark upon the Bunch-Davis (BD) vac¬ 
uum m , which is the minimum energy states, when they 
pop out of vacuum inside the horizon of an inflationary 
background. However, various effects of physics at energy 
scales higher than that of inflation [20] or multi-field ef¬ 
fects [21 can excite these fluctuations to a state other 
than the Bunch-Davies vacuum. In a previous work, we 


2 In models of anisotropic inflation El, generally quadrupole term, 

~ 2 

(k • fi) in momentum space appears which is different from the 
quadrupole term in position space we predict here. 


showed that how by assuming initial conditions other 
than the Bunch-Davies vacuum, one can decrease the 
tensor/scalar ratio in a high energy scale chaotic mod¬ 
els like m 2 (j ) 2 [18] and make it compatible with the latest 
Planck data [lj[22]. We also showed how one can induce 
large amount of running in the scalar spectral index or 
blue tensor spectral index using scale-dependent initial 
condition [23]. It is shown that excited initial state can 
induce larger y -type distortions in comparison with the 
Bunch-Davies vacuum [24] . 

The equation of motion for the gauge-invariant scalar 
perturbations, the Mukhanov-Sasaki variable u{r , y) p 5 j . 


. a 5 (j) 

u = -z (-77 + T 

a cp' 




a 


(5) 


is 




(6) 


Prime denotes derivative with respect to the conformal 
time r and Uk(r) is the Fourier mode of u(r,y). For a 
quasi-de-Sitter background 

a ( r )--A> oo 


where H is the Hubble constant. The most generic solu¬ 
tion to (J6| with ([7]) is 


u k (rj) 



a k HW(k\T\) + 0 k H™(k\T\)] , 


(8) 


where H ^ and H^ 2 are respectively Hankel functions 
of the first and second kind. The terms proportional to 
and /3j£ respectively behave like the positive and neg¬ 
ative frequency modes in infinite past. These Bogoliubov 
coefficients satisfy the Wronskian constraint, 

Kl 2 - lAfel 2 = 1. (9) 


The standard BD vacuum is obtained when ak = 1 and 

Pk = 0. 

For a generic initial state, the energy and pressure den¬ 
sity carried by the fluctuations are of the same order, 
^Pnon-BD ^ ^Pnon-BD, and should remain subdominant 
with respect to the total energy of the inflaton. Also 
their variations with time should not hinder the slow-roll 
inflation. Noting that fyn 0n _ BD ~ Sp' non _ BD ~ Ittfynon-BD 
in the leading slow-roll approximation, this requirement 
is satisfied if 


^Pnon-BD ^ e P0 5 ^Pnon-BD ^ ^P 6 P0 ; 

where e and p are defined as 


( 10 ) 


nr! 

e= 1 - 72 ^ 


(u) 


The strongest of the above two constraints may be writ¬ 
ten in terms of [3k as 


f~ Ki¬ 


rn 
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We will assume that all scales of interest are uniformly 
excited to an initial state with the second Bogolibubov 
coefficient, which is anisotropic in position space within 
the initial inflating patch, 

Pk=P o(£), (13) 


once their physical momenta become smaller than the 
scale of new physics, M [26]. Inevitably, modes which 
remain above this hypersurface momentum do not get 
excited and therefore the left hand side of the integral 
(12) remains finit^] With the choice (13), one does not 
lead to extra /c-dependence in the power spectra and does 
not change the spectral index at the observable scales. As 
mentioned, we also assume that the second Bogoliubov 
coefficient can depend on the position direction through 
the parameter /3q(x). Having 


Vnon-BD ~ |A>(£)| 2 M 3 4 , 5p' non _ BD /K ~ |/? 0 (*)| 2 M 4 , 

(14) 

one obtains the following upper bound on |/?o(£)|, 


\Po(x)\ < 


HMpi 

M 2 


HMpi 

e W 2 


(15) 


As it was discussed in |T 8 ] and will be reviewed briefly 
below, \/3q(x) \ is not necessarily very small. Larger values 
of \/3o(x)\ are compensated with a smaller Hubble param¬ 
eter, iL, for a given model to match the normalization of 
density perturbations with the data. 

The scalar power spectrum defined as, 




_fe^_ 

27r 2 


Uk 2 




fc/SC-K) 


(16) 


turns out to be 


^s = Pbd 7 s > ( 17 ) 


where 

7bd= 8 !^(w) ’ ^ = ( 18 ) 

To study this power spectrum more closely, we note 
that the energy and the power spectra (and also the bi¬ 
spectrum) expressions only depend on relative phase of 
(3 . Hence, they may be parameterized as 


a k = coshx s e* Vs , /?f = sinhx s e . (19) 


3 There is a qualitatively different situation in which at time to , the 
modes with physical momentum smaller than M get pumped to 
an excited state, whereas the larger ones remain in their vacuum. 
These two pictures, even though are qualitatively different, lead 

to the same result quantitatively. For the latter to be relevant 
for the CMB scales, one would expect that inflation did not last 
more than what is needed to solve the problems of Big Bang 
cosmology. 


With this parametrization, ~ sinh 1 /3 0 (£), e 2x s < 
7 s < e 2x s . As it was shown in m, in the regime where 
the deviation from the Bunch-Davies vacuum is large, 
1 , in order to have maximal separation between 
the scale of new physics, M, and the inflationary Hubble 
parameter, iL, one is confined to tp s ~ 7 t/2 . For ra 2 </> 2 , 
for large 1, M ~ 21 H. 

Let us now assume that the horizon patch or the asym¬ 
metric effect of new physics at the energy scale higher 
than the energy scale of inflation is anisotropic and in 
particular it singles out one direction such that the pa¬ 
rameter /3q(x) in the initial state Bogoliubov coefficient 
involves an asymmetric term too. We parameterise this 
asymmetric effect at the new physics hypersurface with 
s as follow^ 


/3o(x) = sinhy s (l + ex.n)e l(p s ; ( 20 ) 


where e < 1. We may interchangeably use cos ^ for 
x • n in the rest of the analysis. From the Wronskian 
constraint one can easily obtain the norm of the first 
Bogoliubov coefficient too. Following the parameteriza¬ 
tion of m the first Bogoliubov coefficient takes the form 

al = [1 + sinhy s (l-b ecos^) 2 ] 1 ^ 2 ^ - ^- 5, . ( 21 ) 

One can easily obtain the factor 7 $ through the relation 


4 In this paper, we assume that the mechanism that excites the 
fluctuations within the horizon patch is position-dependent and 
has picked up a small dipole-dependent correction in addition 
to the usual monopole homogeneous term. This in particular 
is conceivable if one assumes that the horizon patch was bigger 
than a Hubble size and thus the mechanism responsible for the 
excitement of the mode leads to different values in different parts 
of the horizon patch. In the first approximation we assumed that 
the second Bologoliubov coefficient has a small dipole correction 
in addition to the uniform part. The anisotropic vacuum that 
we have hypothesised in our article to justify the hemispherical 
asymmetry would correspond to an anisotropic energy momen¬ 
tum tensor “within” the horizon patch. This should be compared 
with the situation described by Erickcek et al. _9]. where a mode 
larger than our current horizon creates this asymmetry. In that 
scenario, it is not clear why a mono-wavelength superhorizon 
anisotropy in a particular direction should be hypothesised to 
obtain the observed hemispherical asymmetry. The observation 
of the hemispherical asymmetry in the CMB in that scenario 
provides information about the distribution of energy momen¬ 
tum tensor at the scales beyond our horizon. Another variant 
of this scenario which includes all superhorizon modes and non- 
gaussianity on scales larger than the scale of our universe are 
considered Id. tackles the fine-tuning issue better. As we will 
see, what our computations show is that if the energy-momentum 
tensor of the inflaton is anisotropic “within the horizon patch” 
in the beginning of inflation, no matter how long the subsequent 
inflation lasts, the resulted power spectrum is anisotropic. 
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(18). Expanding 7 $ as a function of 6 , one obtains 


7s = 7o + scos Vtr7i + s 2 cos 2 ^72 + • • • 

= cosh 2x s — cos 2(p s sinh 2x s 
+ scos ipx [2 tanh (sinh 2x s — cos 2cp s cosh 2x s )] 

+ e 2 cos 2 'tpx [2 sinh 2 — cos 2 (p s (cosh 2x s + 2 ) 

tanh 3 x s ] +••• 

( 22 ) 


where 7 ^’s are the level i -th order coefficient in e expan¬ 
sion of 7 s and the ellipses stands for the higher order 
of e which are suppressed. The higher order terms are 
not only suppressed by powers of e, but also by powers 
of e~ 2x s where 1 , which make them completely 

negligible in comparison with other terms. In this limit, 
the primordial scalar power spectrum obtains dipole and 
quadrupole directional dependence 


ys = y iso (l + 2A(x • ft) + B(x • h) 2 ) . (23) 


The parameters A and B are respectively 

571 


aennea as 


A = 


B = 


2 70 

£ 2 72 

7o 


(24) 

(25) 


In the large limit, 1, A and B can be expandec 

as 


A = [l — 2e 2x s + 2e 4x s (1 - cot (f 2 s ) + ...] e , 

(26) 

B =[ 1 — 2e -2x s - 2 e _4x s (1 - cot 9 J 2 ) + ...] e 2 , 

(27) 

where ellipses contain terms that are higher order in 
exp(— 2x s ) which may have dependence on ip s too. It 
is interesting that no dependence on ip s appears in the 
leading (and even next-to leading) order of these param¬ 
eters when 1. This is the limit that we we will 

focus on in the rest of our analysis. As shown in [18], 
in this limit maximum separation between the scale of 
new physics and inflationary Hubble parameter, which 
is required to utilize the effective field theory, could be 
obtained. Also in this limit, it turns out that effectively 
cp s ~ In this limit, 


B^A 2 . 


(28) 


One should note that with decreasing , the parameters 
A and B decrease too. Also the relation between these 


two parameters, eq. (28), will not hold any more. 


The Planck data indicates A = 0.072 d= 0.022 on large 
angular scales with the best fit for the anisotropy direc¬ 
tion to be (Z, ft) = (227,-27). This observational con¬ 
straint determines 6 ~ 0.07 =b 0.022 and the h direction 
to be the unit vector along the anisotropy direction. For 
smaller values of x s , the factor 71/70 < 2 , and thus 
one has to increase the required amount of e to account 


for the observed hemispherical asymmetry. The mini¬ 
mum value for 6 , giving rise to the observed hemispheri¬ 
cal asymmetry, is hence ~ 0.07. We conclude that 


5 >0.07 ±0.022. (29) 


In order to induce the asymmetric effect on a finite range 
of scales, one has to assume that e is scale-dependent. 
This in particular can be realized assuming that the 
asymmetric effects become more effective when inflaton 
passes through the scales that left the horizon at very 
large scales. For example if one assumes that the inter¬ 
action between the inflaton and the term that induces the 
asymmetry in the Lagrangian is proportional to \(\> — 0o| n , 
where (j )0 is the value of the inflaton when the scales 
corresponding to our horizon scale left the horizon. As 
the inflaton moves away (j )0 creating asymmetric excited 
quanta becomes more expensive and the effect fades away 
at smaller scales. 

The existence of quadrupole term, B{x • n) 2 with 
B ~ A 2 , besides the dipole term, is one of the predictions 
of the model. With the observed value of A, the coeffi¬ 
cient of the quadrupole term turns out to be quite small, 
0.0025 < B < 0.008. Such would not be discernible from 
the systematics and noise from the current CMB data. 
Higher order multipole coefficients are also present in the 
primordial spectrum, but suppressed by the correspond¬ 
ing power of e. This is a typical feature also among var¬ 
ious previously proposed power asymmetry -generating 
mechanisms, that the predicted couplings of the CMB 
angular modes falls of with their multipole separation. 
However, for example, in the non-Gaussianity -generated 
cases of Ref. [12], the suppression seems to be square- 
root rather than linear e^/ 2 , less sharp than in our case. 
Higher multipoles can thus be used as a cross-check to 
distinguish between different explanations for the origin 
of the dipole asymmetry. 

In the next section we will investigate another signa¬ 
ture of the model in the bispectrum, which might be 
easier to detect. 


III. BISPECTRUM 


Let us calculate the three-point function for the above 
direction-dependent excited states to see how they mod¬ 
ify the bispectrum. One can calculate the Wightman 
function for the solution, Uk(r), 


Gk (r, t ') 


-H -2 Mfe(r) u* k (r') 
02 a (r) a(r') 


(30) 







The three-point function could be derived from the 
Wightman function through the following integral [29] : 


5 


(E^) M ^ H l° o dT 

T 2 (a(r)d T G^ (0, T))(a(r)d T G^ 2 (0, T))(a(r)d T G^ 3 (0, r)) 
k 3 

^permutations + c.c., (31) 

where To is the moment at which the physical momentum 
becomes equal to the physical cutoff, To = The 

Wightman function in the integrand is 



drGfc (0 ,t)) = 


H 3 2(a t -fe)K e jfcT + /?fe e~ ifcr ) 


2 4> 2 k a{r) 

The bispectrum takes the form 


(32) 


2# 6 E 


^kSk 2 ^k 3 ) — 


<i>j 


i= 1 


j> 2 M 2 P n(2fcf) 

i=1 


^ 1 - cos(fc t r 0 ) ^, sin(fc t r 0 ) ^ 1 - cos(fcjT 0 ) 

h-L Lg 3 


3 =1 


kn 


-E- 

i=i 


sin(fc :/ To) 
? kn 


(33) 


FIG. 1. The amplitude of non-gaussianity depends on which 
corner of the sky the mode is located. In a scale-dependent 
asymmetry, it also depends on whether the asymmetric effects 
are effective at the scale of interest. 


function one obtains: 


3 H s e 

* — 16 ( 27 rfs^ki)—, -^-e, 


i=1 0 2 n(2 kf) 

2=1 




(35) 


where 


where k t = &i + &2 + and kj = k t — 2kj. Terms propor¬ 
tional to ^ and are respectively the ones that can lead 
to enhancement in the local configuration, k\ ~ /c 2 &3 

[30] , or flattened (folded) configuration, k\ + — &3 [31] . 

Coefficients, 3 ^, ^ and are as follows 



: n (“fci -^j(ii a 

1 + 


+ c.c. 

£8 = *II(a fei - 3ki)(— II a 

1 + 

n«.) 

+ c.c. 


-<fc.)(§ t IK + 

Ri u . 

a kj 

K 

n^) 

+ c.c. 

= * R( Q; fti 

kj 

at. 

Kj 

n^) 

+ c.c. 





(34) 


C = Re 




n«,ff+fi' 


,v v«. 

3 

= Re II(^ - &«)« + + <4^) 


. 2=1 


(36) 


One can calculate the /nl parameter using the defini¬ 
tion, 


/nl = - 


5 s(c kl Ck 2 Ck 3 ) 

6 J^iCkiCk^iCkjCkj) 5 

2 >3 


(37) 


The enhancement of the flattened configuration is how¬ 
ever lost in slow-roll inflation after the projection of the 
bispectrum shape on the 2-dimensional CMB surface |32j • 
Besides for the large deviations from the Bunch-Davies 
vacuum where, 1 and (j) 7t/2, the enhancement 

factor is exactly equal to zero.Thus we focus on the local 
configuration enhancement. Noting that k\ ~ & 3 , 

the enhancement for the local configuration three-point 


to be 


/nl — - 


20e ki 


3 k 3 7s(&3)(7s(&i) + ls{k 2 )) ' 
Expanding the £/nl in terms of 5 up to second order, 

/nl - /nl + /nl £ + /niE - 


(38) 


(39) 



















6 


I k 3 l«l k x W k 2 M k l=2500 I 

I^W^iol 



I & 3 l«l ^ W k 2 M k l=2500 I 
I k 3 l«l & /=10 I 



FIG. 2. The amplitude of local non-gaussianity reaches its 
minimum departure from the mean value of non-gaussianity, 
when the considered modes are around the preferred pole of 
the sky and \k 3 \ ~ |fe«io| \ki\ ~ \k 2 \ ~ \ki ~ 2500 1 • This is 
because for e{k\) — £{k 2 ) = 0 and £(k 3 ) = £ ^ 0. 


FIG. 3. The amplitude of local non-gaussianity reaches its 
maximum departure from the mean value of non-gaussianity, 
when the considered modes are around the antipode of the 
preferred pole and \k 3 \ ^ |fe~i 0 | <C \ki \ « \k 2 \ « |fe~ 2 50 o|- 


in the limit that Xs 1 an d th us Ps — 77we have 

< 4 °> 

/nl - y ^ [cos(V'* 1 ) + cos{ipz 2 ) - 2 cos(V>£ 3 )], (41) 

/nl - -y ^ [cos(V’a 1 ) 2 + cos(^j. 2 ) 2 - 6 cos(^ S3 ) 2 

4 cos( / ip £ 1 ) cos('ip £ 2 ) + 4cos(V’x 1 ) cos(V>x 3 ) 

+ 4 cos^) cos( , 0x 3 )] • (42) 

Hence the amplitude of the bispectrum depends on the 
angles that three position vectors make with the preferred 
direction, please see fig0 

The first term, which gives the dominant con- 

tribution to the bispectrum, is independent of the these 
angles though. For an inflationary model with e « 0.01, 
the 

4l-4.17, (43) 

where we have taken the largest scale at which the cosmic 
variance is negligible to be corresponding to / m i n ~ 10 
and the smallest one to be the largest l probed by the 
Planck experiment, / max — 2500 [^] This is within the 2a 
allowed region of local nongaussianity from the Planck 
2015 experiment |33j. At higher orders, the excited asym¬ 
metric initial condition induces directional dependence 
to the bispectrum at the first order correction. It is easy 
to verify that if X\ = x^ = x% and asymmetry is scale- 
independent, the directional dependence vanishes at both 


5 There is a minor correction to to the value of /nl as quoted in 
m due to a missing factor of ~ 10 in the previous analysis. 


first and second order in 6. That is if the chosen momenta 
are all at the same corner of the sky, there is no modu¬ 
lation on top of the mean value (43). On the other hand 


if the three modes are at different corners of the CMB 
sky, even if the asymmetry is scale-independent, one will 
see modulation on top of the mean non-gaussianity value 
(43). The maximum of the modulation occurs when the 


short wavelength modes, k\ and are at the preferred 
pole, cos('0^ 1 ) = cos(Vt£ 2 ) = 1 and the long wavelength 
mode is at its antipode. For e = 0.07, the maximum 
non-gaussianity obtained at the scales is 


/nl* - /m( 1 + 4e + 10e 2 ) « 5.54. 


(44) 


On the other hand, the minimum shift from the Bunch- 
Davies value result, corresponding to minimum value for 
non-gaussianity, would occur for when the short wave¬ 
length modes are the antipode pole and the long wave¬ 
length mode is at the preferred pole. For this configura¬ 
tion 


Mt - f { NL( 1 - 4e + Kk 2 ) « 3.20. (45) 


Since the hemispherical asymmetry parameter is scale 
dependent and fades aways at l > 600, one should recon¬ 
sider the above result. In fact the first order correction to 
the mean value of non-gaussianity, which is the dominant 
modulation term, could be written as 


A A 1 ) — 

5 6 k\ 

— — [e(fci) cos(f/>s 1 ) + s(k 2 ) cos(ipx 2 ) ~ 2er(fe 3 ) cos(^ 3 )] , 
o 

(46) 


maximum modification from the mean value (43), occurs 
for the configuration where k\ ~ k 2 is the wavenumber 
corresponding to l ~ 2500 and k% corresponds to l = 10. 
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In this case, e(fci) = e(k 2 ) = 0 and e(ks) = 5 . Now if the 
large wavelength mode, k% is around the preferred pole, 
cos (/0x 3 ) = 1, we will have a decrement from the mean 
value of non-gaussianity 

/Si n ^/t°hl-2e + 3 £ 2 3 4 5 6 7 )« 3.64. (47) 

It does not matter at which corner of the sky, k\ and &2 
modes are located. We can assume that they are centered 
around the preferred pole too. 

On the other hand, If k% is around the antipode of the 
preferred pole, there will be an enhancement in the mean 
value of local nongaussianity 

/NL ax - Cl (1 + 2e + 3e 2 ) « 4.81. (48) 

The difference between these two values of non¬ 
gaussianity at two poles, £/nl ~ 1.17, can be used to 
distinguish this scenario from the competing proposals 
that try to explain the observed hemispherical asymme¬ 
try. 

We recall non-gaussianity can increase the probabil¬ 
ity of a power asymmetry. The impact of higher order 
nongaussinity has been also considered Q2CE3]. a high 
#7VL, even if occurring in an isotropic vacuum, could in¬ 
duce anisotropic ” power asymmetry” on the /nl- It is 
an open question at the moment, whether a gNL, while 
staying within the observational bounds, could induce 
anisotropies in the /a tl at the order we are predicting. 


IV. CONCLUSION 

Hemispherical asymmetry, is one of the persistent 
forms of isotropy violation that has been observed in both 
WMAP and Planck data. The hemispherical asymmetry 
observed is scale-dependent and vanishes for l > 600. 


In this paper we suggested a scenario that accounts for 
the observed hemispherical asymmetry using asymmet¬ 
ric initial condition. We noticed that there are infinite 
higher multipole corrections to the power spectrum. In 
the limit where the scale of new physics is maximally sep¬ 
arated from the Hubble parameter, only the dipole and 
quadrupole terms survive. Requiring that the amount 
of asymmetry is as large as the observed value, we ob¬ 
tained a small but finite amplitude of quadrupole correc¬ 
tion to the power spectrum. This is one way to distin¬ 
guish this scenario from other scenarios that try to ex¬ 
plain the origin of hemispherical asymmetry. We showed 
that the model exhibits unique signatures in the bispec¬ 
trum. Due to the excited initial states, the local config¬ 
uration gets enhanced around a mean value which for 
a large field inflationary model, e ~ 0.01, turns out 
to be around 4.17. There will be modulations on top 
of this mean value which is dependent upon the an¬ 
gle that the patches that contain the mode make with 
the preferred direction. The amount of variation mini¬ 
mizes for the configuration in which the short and long 
wavelengths modes are ^around the preferred pole and 
\h\ « |fcj«io| < |£i| ~ 1*21 ~ |£ j « 250 o | with /^ x « 4.81. 
The maximum occurs when these modes are at the an¬ 
tipode of the preferred pole, /^l 1 ~ 3.64. The difference 
of non-gaussianity between these two configurations is as 
large as ~ 1.17 which would be a definite indication of 
the model. 
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